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Abstract 

The electronic structure of the cubic and (high pressure) hexagonal phases of EU2O3 have been 
investigated by mean of full potential linearized augmented plane wave calculations, within the 
LDA+C7 method. For the cubic phase, the comparison between ferromagnetic and antiferromag- 
netic calculations shows that the exchange interaction is very weak and is therefore expected to 
have a negligible effect on the magnetic susceptibility. This is consistent with the experimental 
behavior of the susceptibility of solid solutions of EU2O3 into A2O3 (A=Y, Lu, Sc). The calcu- 
lations performed for the high pressure hexagonal phase, on the other hand, show that there is 
an antiferromagnetic exchange interaction between nearest neighbor Eu ions, which should have a 
sizeable effect on the susceptibility. Our results allow us to discuss the existing theories. 
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I. INTRODUCTION 



Rare earth compounds present a large variety of interesting magnetic behaviors, due to 
their partially filled /—shells. Trivalent europium compounds, in particular, offer an unique 
chance of observing exchange coupling between " non- magnetic" ions. The ground state of 
Eu 3+ (4/ 6 ) ions is 7 Fq, with a total angular moment J = resulting from both L = S = 3 (in 
atomic units), however, spin- dependent exchange effects are present. Furthermore, Eu ions 
show a substantial admixture of higher energy J = 1 states, which contributes significantly 
to their susceptibility (%). 

EU2O3 is the prototypical compound of this family, its magnetic susceptibility has been 
experimentally investigated in detail^-. The first attempt to explain its magnetic behaviour 
was made by Huang and Van Vleck-, who showed that the susceptibility of EU2O3 is larger 
than the corresponding susceptibility of the free Eu 3+ ion, because the energy levels of 
the excited 7 Fi states are modified by the crystal field, and pointed out the dominant 
contribution of the Van Vleck component of the susceptibility, Xvv- An explicit calculation 
of the susceptibility from energy levels obtained by optical spectroscopy measurements^ 
was performed in this work. However, the resulting \ turned out to be smaller than the 
experimental value, and so the remainder was attributed entirely to the exchange coupling 
among Eu 3+ ions^. 

This viewpoint has been challenged on the basis of measurements of x i n solid solutions 
of EU2O3 into A2O3 (A=Y,Lu,Sc)^ 13 -^, where Eu atoms are diluted. As the number of Eu 
nearest neighbors is reduced, the total interatomic exchange interaction should decrease and 
become negligigle at small concentration of Eu, and the total susceptibility should decrease 
along with it. The experimental susceptibility (per mole of Eu), however, does not decrease 
in the full range of investigation (i.e. up to 10 % of EU2O3 into A20jr> 3 -^), suggesting a 
negligible role of the exchange interaction. This decrease of the exchange component on 
dilution might be due to compensation caused by change in the crystal field splitting of Eu 
sites in A2O3 in comparison with EU2O3; this explanation, however, is not supported by 
optical measurements of the energy levels^*^. Yet another explanation has been proposed 
based on the distribution of Eu atoms in the two available sites (symmetry S6 and C2) of 
the cubic bixbyite structure of these oxides^. The decrease of the exchange component 
on dilution might be compensated by a preferential occupation of the site in the solid 
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solution; the Van Vleck susceptibility of this site is larger because the Fi levels are lower 
in energy in comparison with the C2 site. This preferential occupation, however, was not 
found in X-ray diffraction^, and Mossbauer spectroscopic studies^ 1 ^. 

A further criticism to viewpoint that the exchange coupling guides the physics of su- 
ceptibility also comes from the calculations of the Van Vleck susceptibility performed by 
Caro et al.— . These authors determine the crystalline field parameters starting from the 
experimental values of the energies of the excited states, and the matrix elements of the 
Van Vleck susceptibility were calculated using this potential and the atomic functions. The 
susceptibility calculated accordingly agrees well with the experimental one^, validating a 
picture in which the exchange coupling should be negligible. We should mention here, how- 
ever, that the optical data used by Caro et al. have been subsequently corrected by more 
recent measurements^. 

Therefore, the question of whether the excess of susceptibility may be entirely attributed 
to the exchange coupling is still open. Absence of long range magnetic order does not allow 
direct evaluation of the exchange coupling constants starting from experimental data and so 
a theoretical determination of these parameters by means of ab-initio calculations is highly 
desirable for understanding the behaviour of susceptibility in the material. Remarkably, to 
the best of our knowledge, there is no experimental evidence in the literature of exchange 
coupling between J = ions making the "excess" of susceptibility of EU2O3 a novel puzzle 
where ab-initio calculations have the potential of providing interesting evidence towards the 
solution of this puzzle. 

From the theoretical point of view, rare earth compounds represent a challenge for modern 
electronic structure calculations. Their multiplet structure cannot be explained on the basis 
of a single Slater determinant. The traditional density functional theory (DFT) methods 
with local/semi- local approximations to the exchnage correlation functionals fail to describe 
their correlated nature and result into a qualitatively wrong picture with flat / bands ac- 
cumulated around the Fermi level (Ep). To overcome these problems the LDA+C7, and 
the self-interaction corrected local density approximation (SIC-LDA) methods have been 
widely used in the past. Focusing entirely Eu 3+ compounds, LDA+U calculations were per- 
formed by Johannes and Pickett^ on EuN and EuP, and by Deniszczyket al.— on EuF 3 and 
EUC02X2 (X=Si,Ge). SIC-LDA calculations were performed on several rare earth oxides by 
Petit et alM*H. All these works show that the correct physics in these materials can be 
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treated within the LDA+U (SIC-LDA) method. Hence in the present work we compute the 
electronic and magnetic structure of EU2O3 using the LDA+U approach. Eu being one of 
the heavy rare earths, the use of an all-electron method is almost essential- in the present 
work we emply full-potential linearized augmented plane wave (FLAPW) method^ imple- 
mented within the Elk code^. From the corresponding results we obtain very small values 
for the exchange coupling constants in cubic EU2O3, which imply a negligible contribution 
of exchange to the magnetic susceptibility. 



II. COMPUTATIONAL DETAILS 



For the cubic structure a k-point mesh of 2 x 2 x 2 is used for the Brillouin zone (BZ) 
integration; given the large dimension of this system (40 atoms per unit cell) and the insu- 
lating character of the compound, this choice is a reasonable compromise between accuracy 
and computational load. For the (smaller unit cell) hexagonal structure, a mesh of 4 x 4 x 4 
is used. Spin-orbit coupling has been included in all self-consistent calculations. We have 
used the fully- localised-limit of the LDA+U method^. The values of U for Eu /-states and 
J are chosen to be 7 eV and 0.75 eV respectively^. 



III. CRYSTAL AND MAGNETIC STRUCTURE 

EU2O3 crystallizes in the cubic bixbyite structure (shown in Fig. [[]) with space group 
(ia3); it corresponds to a bee lattice with a primitive cell of 40 atoms^. Eu atoms occupy 
the two inequivalent Wyckoff sites 8a (site symmetry Sq) and 24d (site symmetry C2), with 
a distorted octahedral coordination of O atoms (in the 48e sites). The experimental lattice 
constant for EU2O3 is a = 10.859 A—. The bixbyite structure is derived from a defective 
cubic fluorite structure, which corresponds to a simple cubic lattice with 10 atoms and half 
value of the lattice constant. Eu atoms occupy the positions of a fee lattice; the defective 
structure is obtained by removing two O atoms at 1/4 and 3/4 of the body diagonal. Due 
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to the large number of atoms in the unit cell in the bixbyite structure, Ref. 
smaller, defective fluorite unit cell. In order to compare our results with this work, we have 
also performed calculations within the defective fluorite unit cell. 

In order to evaluate the magnetic exchange coupling constant, the band calculations 
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FIG. 1: (Color online) (a): Spin structure of the E^Os determined using neutron diffraction^ 3 - 
method, (b): Full structure of cubic EU2O3. (c) O coordination around Eu 56 and(d) O coordina- 
tion around Eu C2 sites. 

need to be performed with different configurations of the Eu spins. Experimentally, EU2O3 
does not have a long range magnetic order; thus as a starting point we have used the 
spin configuration of isostructural Er 2 03 compound (see Fig. [1]). This spin configuration 
also turns out to be a stable for EU2O3 in our calculations. However, in this Er 2 03-like 
configuration the sum of the scalar products of the spin of one ion with its nearest neighbors 
(NN) spins is always zero, which does not allow for a determination of the NN coupling 
within a Heisenberg model. In order to obtain a positive sum of scalar products one half 
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FIG. 2: (Color online) Left panel: high pressure hexagonal structure of EU2O3. Right panel: 
antiferromagnetic configuration of the hexagonal phase. 

of the spins of the C2 sites have been reversed relative to the E^Os-like configuration such 
that all the spins point along the positive direction of the cartesian x, y or z axes. In order 
to compute the exchange couplings we have also studied other ordered phases- to obtain 
a ferromagnetic (FM) like state where the direction of the Sq sites is chosen to maximize 
the number of positive components, which may be written as (1,1,1), (1,1,-1), (-1,1,1) and 
(1,-1,1) in the appropriate units; we obtain in this way a positive sum of scalar products for 
a Sq or C2 europium with its twelve NN. In order to obtain an antiferromagnetic (AFM) 
like state, the Sq spin direction is reversed; the sum of products is then negative. 

EU2O3 undergoes a structural transition under pressure from the cubic to the hexagonal 
phase (space group P3ml)^ (see Fig. [2]), which has 5 atoms per primitive cell (Eu atoms 
on the Wyckoff positions 2d and O atoms on the 2d and la positions)^. The transition 
begins at 5 GPa and completes at 13 GPa. At this pressure, the lattice parameters are 
a = 3.738 A and c = 5.632 A—. For this structure, the band calculations are performed 
with FM as well as AFM spin configuration; in the AFM configuration (shown in Fig. [2]), 
the two cations of the primitive cell have antiparallel spin orientation. 
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IV. RESULTS AND DISCUSSION 



Before calculating the electronic structure a full structural relaxation is performed and 
the results are compared with the available experimental data. Due to the large size of 
the bixbyte unit cell, these calculations for the cubic phase have been performed using the 
defective ffuorite cell with a FM configuration. The volume dependence of the energy was 
fitted to the third order Birch- Murnaghan equation of stated. We obtain an equilibrium 
lattice constant ao = 5.311 A, 2 % smaller than the experimental value (consistent with the 
typical error in LDA or LDA+[7 calculations). The calculated bulk modulus B = 140 GPa, 
is also in good agreement with the experimental value of B Q exp = 145 ± 2 GPa^ (measured 
in the bixbyte phase). 

As mentioned earlier, the first problem in the study of this system is the choice of spin 
configuration, since experiments do not report any magnetic order in this J = system. We 
therefore assume as a starting (reference) configuration the experimental magnetic structure 
of isostructural Er 2 03. In this configuration, spin moments are not collinear, and the two 
crystallographic sites have totally different directions^: while the Sq moments are directed 
along the diagonal axes, C2 sites direct their spin along the Cartesian axes (see Fig. [T]). 

The band structure of cubic EU2O3 with the E^Os-like spin configuration is shown in 
Fig.[3j The red squares and blue circles in Fig. [3] represent the importance of the contribution 
from Eu 4/ states, for sites S$ and C2 respectively. The general features of EU2O3 bands 
can be readily explained as follows: the O p states are responsible for the bands located 
in the energy region between w —3 eV and the Fermi level (which is the zero of energy 
in all our plots). The Eu 3+ ions are in a 4 /^configuration, which leaves an empty 4/ 
orbital per atom in the majority spin channel. Minority 4/ states, on the other hand, are 
completely empty and form the group of bands from +4 to +8 eV. Occupied majority spin 
4/ states form the group of 24 bands sitting from -5.5 to -3 eV and their splitting into four 
separated groups is similar to the situation in Eu pnictides^. This splitting in the case of 
Eu pnictides was interpreted as an effect of the intra-atomic anisotropic exchange- a test 
calculations based on LDA resulted into a single narrow /-manifold indicating that splitting 
cannot originate from crystal field effects, making the complex exchange effect in this open 
shell system the most likely explaination. The large distance between (equal spin) filled 
and empty counterparts derives from the large U value for localized 4/ states. The empty 
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FIG. 3: (Color online) Band structures of cubic EU2O3 along the symmetry directions of the 
Brillouin zone. Eu 4/ character is denoted by red squares and blue circles for the Sq and C2 sites 
respectively. 

majority spin states lie around 2.3 eV above Ep. This electronic structure therefore reflects 
the electronic configuration of the Eu 3+ ions. All the 4/ bands have a negligible dispersion 
0.2 eV), consistent with the localized nature of these orbitals. 

For better understanding of the electronic states presented in Figs. H] are the densities 
of states (DOS) of cubic EU2O3 in the Er 2 03-like spin configuration. Figs. H] confirms the 
picture given above. In particular, the states from -3 eV to Ep derive from O states, and 
hybridize only weakly with Eu states. Besides 4/ states, Eu's most prominent contribution 
to the electronic structure is the bd states, mostly located in the conduction bands region 
above 4 eV. The narrow peaks from -5.5 to -3 eV and the structures above 4 eV are derived 
from filled majority and empty minority Eu 4/ states respectively. The peaks in the region 
w 2 — 2.5 eV above Ep are predominantly Eu states (in particular, empty majority spin 4/ 
states), with a very small contribution from oxygen orbitals. The marginal Eu 4/- O 2p 
hybridization indicates towards a small value of the NN exchange coupling. 

Importance of the crystal field (CF) effects in this compound has been a subject of 
discussion in past^. Optical measurements have been interpreted assuming that the £6 and 
C 2 sites have a very different CF splitting. Our results support this interpretation (see Fig. 
H]) and show that the / states behave differently for the two sites- the lower symmetry of the 
C2 site results in a larger subdivision of the 4/ peaks with occupied / states shifted 0.5 eV 
lower in energy compared to the occupied /-states with S$ site symmetry. The unoccupied 
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FIG. 4: (Color online) Orbital projected DOS of cubic EU2O3 per primitive cell in the Er203-like 
configuration. The Eul and Eu2 labels indicate the Sq and C2 sites respectively. 

/-states show a similar behaviour with Ci site symmetry states moved higher in energy with 
respect to Sq site projected /-states. This leads to C 2 site having a much larger energy gap 
between occupied and empty majority spin states as compared to Sq site, which in turn 
implies different density matrices for the two sites. This difference can be related to the 
different expectation values of orbital momentum (to be discussed later). 

In order to understand the exchange interaction in EU2O3 we studied two further spin 
arrangements, refered to as FM and AFM configurations (we refer to Sect. Ill II for their 
description). In Fig. [5] we compare the DOS for EU2O3 in the FM and AFM configurations. 
It is clear that spin ordering affects the Eu states marginally and the corresponding DOS 
are almost identical. These results have consequences on the exchange parameters as the 
total energy difference between these two configurations is very small. 

We have also calculated the expectation values of the quantum operators L, S and J for 
the two different sites. The quantization axes are different for the two sites: in the case of C2 
it points along the Cartesian axes, and in the Sq site along the main cube diagonal. In units 
of h we obtain for the Sq site S a = 2.77, L a = —1.56 and J a = 1.21. Our value of S a differs 
from the expected value S a — 3 for the free ion probably because of the use of muffin-tin 
sphere for computing integrals, and also because of an incomplete spin polarization. In the 
case of the C2 site, on the other hand, we get very different results: the components along 
the corresponding local quantization axes a' are S a > = 2.41, L a i = —2.34 and J a > = 0.07. In 
other words, the J = ground state of Eu seems to be reproduced to a much larger extent. 
Again (consistent with experiments) we find very different behaviors for the two sites. We 
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FIG. 5: (Color online) Total DOS of cubic EU2O3 in the Er203-like configuration, and in the two 
(FM and AFM) configurations described in the text. 

should mention that our small value of L a for the site is similar to the value L a = —1.5 
obtained in EuN by Johannes and Pickettr^, who ascribed this result to an overquenching 
of angular momentum caused by an overestimated crystal field effect within DFT. 

At this point it is interesting to compare the electronic structures of EU2O3 and EuN— 
(which crystallizes in the rocksalt structure and correspond to an Eu 3+ configuration). The 
general features of the band structures are quite similar, a part from obvious differences 
associated with the different ligand (the center of gravity of O 2p states lies slightly deeper 
in energy). A qualitative difference exists in the properties of EU2O3 and EuN: EuN is 
metallic in nature and this metallic character arises from partially occupied dispersed Eu 5d 
bands, crossing the empty majority 4/ bands. On the other hand in case of EU2O3 similar 
dispersed bands with Eu 5d contribution exist but they also have a relevant interstitial 
character and lie above the empty majority 4/ bands; as a consequence, EU2O3 turns out 
to be an insulator. It is interesting to notice that both in our calculations for the simplified 
structure suggested by Petit et al.— and in their SIC calculations these bands cross Ep 
and in contrast with experiments gives a spurious metallic character to EU2O3. The origin 
of the difference in character between the bixbyte and the simplified structures lies in the 
relaxation around the vacancies, which is not allowed by symmetry in the latter. In fact, 
an inspection of the charge corresponding to these dispersed bands in the model structure 
shows that these bands derive from quantum states localized in the O vacancy sites. 

Under pressure EU2O3 undergoes a structural phase transition to a hexagonal phase (see 
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FIG. 6: (Color online) Band structures of hexagonal EU2O3 in the primitive cell. Left panel: 
ferromagnetic (FM) state. Right panel: antiferromagnetic (AFM) state. Eu 4/ character is denoted 
by red circles. 

Fig. [2]). Due to the very large computational load involved, the calculation of the transition 
pressure is out of the scope of the present investigation. However, we minimized the total 
energy of EU2O3 in the hexagonal phase. The fitted energy minimum is at Vo = 67.62 A 3 , 
and Bq = 180 GPa. The calculated bulk modulus compares well with the experimental 
value -Bo.exp = 151 ± 6 GPa 24 . As for the equilibrium volume, the experimental hexagonal 
volume is V exp = 68.15 A 3 when the transition is completed (at 13.1 GPa) 24 . Using the 
theoretical bulk modulus, we arrive at an equilibrium volume at 13.1 GPa equal to 62.7 A 3 , 
which underestimates the experimental lattice constant by about 2.8 % (which is typical of 
LDA based functionals). 

In order to understand how pressure modifies the electronic structure of this compound, 
we have studied its properties in the hexagonal phase, at the lattice parameters correspond- 
ing to P = 13.1 GPa. The band structure of the hexagonal EU2O3 is shown in Fig. [6] for 
both FM and AFM configurations. The general features of the bands are similar to the 
cubic phase; hexagonal EU2O3 is a semiconductor, with a gap of 2.3 eV between valence 
and unoccupied 4/ band, and a gap of 4.1 eV between valence and conduction band. This 
semiconducting behaviour is in agreement with the experimental optical and transport prop- 
erties of EU2O3. In contrast to the cubic phase, in the hexagonal phase the dispersed band 
with the minimum around 3.5 eV at the T point are not present. This difference may be 
related to the fact that the hexagonal structure has no vacancies in the atomic O positions 
ruling out the possibility of having bands with a large interstitial character. 

The total DOS of cubic and hexagonal EU2O3 are presented in Fig. [71 While the structure 



11 



14 r 
12 - 

3 io - 



5 8 - 



8 4- 





FIG. 7: (Color online) Total DOS of cubic and hexagonal EU2O3 in the FM state. 



of levels is similar, the O 2p states have larger bandwdith and overlap with the occupied Eu 
4/ manifold. This increase of the O 2p bandwitdh is clearly a consequence of the compression 
O atoms experience under pressure. Also, the exchange-related splitting of the 4/ levels is 
different in the cubic and hexagonal phases. 

Some results of our calculations on cubic EU2O3 may be compared with experimental 
values. The gap between the valence and the conduction band (AE V _ >C ) has been obtained 
by optical measurements; the calculated value AE V _ >C = 3.50 eV is close to the experimental 
one AE opt = 4.3 ±0.3 eV— . The gap between the valence and the empty majority spin Eu 
/-band has been determined by the experimental curve of conductivity vs. temperature; the 
calculated value AE v _ > f = 2.18 eV is in reasonable agreement with the experimental value 
of AE cond = 1.84 eV-21. 

Our calculations allow us to calculate the electron density on the nucleus, which is relevant 
in the context of the Mossbauer spectroscopy of 151 Eu in cubic EU2O3. The electron density 
on the nucleus may be evaluated in terms of the difference relative to tetragonal EUF3, which 
is the reference compound. Using the experimental isomer shift 5 = 1.03 ±0.01 mm/s^ and 



the differential nuclear radius of Ref. |28|, we find an experimental density difference to be 
Ap(0) = 3.18 ± 0.03 ag 3 ; the calculated value for the dominating C2 site is Ap(0) = 3.17 
aQ 3 which is in very good agreement with the experimental value. 
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V. EXCHANGE INTERACTIONS 



The exchange interaction in ions with L ^ is characterized by the dependence of the 
exchange integral on orbital orientation. According to Van Vleck and Huang22 this effect 
gives rise to an "anisotropic exchange", resulting from the dependence of orbital charge 
density on the direction. The coupling between the ions i and j may be described by the 
exchange potential^ 

k, = -EEEE E w^iOW^(2 + 2 *- s j) C 1 ) 

i=l j = l fj,,fi' u=-fi v' =—p! 



where Yn V are the tesseral harmonic operators equivalent to Ref. 



29; rij is the number of 



electrons in incomplete si 
As demonstrated in Ref. 



rell of the ion i, Sj is the electron spin and /i = // = 6 for Eu 



3H 



29l . the exchange coupling for Eu + ions in their ground state and 



in cubic compounds, may be written in the standard form 

V ex = — 2a,gjj-Sj ■ Sj (2) 

where Sj is the spin of the ion % and a^jj are the effective exchange constants. In other 
words, we can use the standard form of the isotropic coupling, to deal with an anisotropic 
exchange. Therefore, it is possible to determine the a^jl constants by calculating the energies 
of different spin configurations. These constants allow to determine the effect of the exchange 
coupling on the magnetic susceptibility. 

Our results show that there is no significant total energy difference between the FM-like 
and the AFM-like phases (for their definition see section UTTj) . with an energy convergence 
parameter e~3x 10~ 4 eV and a total energy E t ~ 5 x 10 6 eV. The maximum possible 
strenght of the exchange coupling consistent with these results corresponds to an energy 
difference per primitive cell AE t ~ 2e which leads to a difference per Eu atom of AE ~ 0.04 
meV. In this work we consider only the exchange coupling of the Eu 3+ ion with its twelve 
NN, since it is expected to be the largest. In fact, in the bixbyite structure, NN exchange 
is mediated by the oxygen atom, which is not the case for next NN Eu atoms. With this 
assumption, the effective exchange constants refer only to the coupling of NN Eu ions. We 
also assume that all the NN pairs have the same a^jj = a e ff- Starting from the energy 
difference AE ~ 0.04 meV and the spin values given by our calculations, we obtain that the 
upper bound of the effective exchange constant to be a e // ~ 0.002 meV. 
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Due to lack of long range magnetic order—, the value of a e // cannot be verified by 
direct comparison with experimental magnetic data. However, this interaction may give a 
contribution to the magnetic susceptibility. According to Ref. |4|, the magnetic susceptibility 
X may be written as 

X = Xdia + Xp = Xdia + XVV + Xex (3) 

where the paramagnetic susceptibility Xp, obtained by subtraction of the diamagnetic core 
component Xdia, is the sum of the Van Vleck contribution xvv and of the exchange contri- 
bution XeA At T = K, 

8iV/i| 3 1 , A . 
Xvv = E p- ( 4 ) 

128N^ B A eff 

Xex ~ Ei(Ei - 164*/) [b) 
where N is the number of atoms, {ib is the Bohr's magneton, K B is the Boltzann's constant, 
are the energies of the triplet state 7 Fi centered around the energy Ei, and 

Aeff = E a { :}} = 12a eff (6) 
i=i 

However, the contribution of the Sq and C2 sites must be evaluated separately, because the 
energies of the 7 F! state are different. Since the multiplet energies are not accessible from our 
one-particle calculations, we calculate the Van Vleck contribution by using the experimental 
optical energies for the two sites given by Ref. [(3; we get Xv^ = 7-72 X cm 3 /mol 
of Eu (in CGS units). If we subtract this contribution from the experimental value of the 
paramagnetic susceptibility x P exp ^ = 9.1 x 1CT 3 cm 3 /mol of Eu (in CGS units) of Ref. [j, we 
obtain the value of the excess susceptibility Ax = X P ~ Xvv = 1-38 x 10~ 3 cm 3 /mol of Eu. 

In accordance with Van Vleck and Huang^22 if we assume that this excess susceptibility 
is due to the exchange contribution only (Ax = Xex) and that A e ff is equal for both sites, 
we can estimate the experimental value of A e ff = 0.41 meV. By using the maximum value 
of a e ff obtained by our band calculations, we obtain a theoretical upper bound for the 
constant to be A e ff ~ 0.02 meV. Our results, therefore, lead to the conclusion that in 
cubic EU2O3 the difference between the experimental paramagnetic susceptibility and the 
Van Vleck contribution cannot be due to the contribution of the exchange interaction, in 
contrast with conclusions of Van Vleck and Huang^22. It is consistent, on the other hand, 
with the fact that the experimental susceptibility (per mole of Eu) does not decrease in 
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the solid solutions of E112O3 into A2O3. The point of view of Van Vleck and Huang has 
also been criticized before by Caro et al.—, who performed a calculation of the Van Vleck 
contribution xvv including the matrix elements among atomic wavefunctions of the Eu 4/ 
states. They obtained a resulting Xp m good agreement with experiment, without invoking 
any contribution from the exchange interaction. Our calculations are also consistent with 
these results. 

The investigation of the exchange interaction has also been performed for the hexagonal 
phase of EU2O3 under pressure. In this case, Eq. [2] is not applicable, because the structure of 
the compound is not cubic; therefore, it is not possible to determine the exchange constants 
by total energy differences. In contrast to the cubic phase in the case of the hexagonal 
structure the application of a weak magnetic field along the c axis gives collinear spins. In 
the AFM configuration (described in Section IIIII) the cation spins have parallel orientation 
in planes perpendicular to the magnetic field direction (z axis), with antiparallel orientation 
between planes. The Eu atom has three neighbors at 3.50 A in the lower plane, three at 3.59 
A in the upper plane, and four neighbors at a distance of 3.74 A in the same plane^ 1 ^. The 
energy difference per Eu atom is Eafm — Efm = —3.63 meV. Therefore in the hexagonal 
structure we obtain a small but sizable AFM exchange interaction between the Eu ion and 
its six NN. Interestingly the energy difference is of the same order of magnitude as that 
found by Johannes and Pickett for EuN and EuP— . 

The exchange interaction is about 100 times stronger in the hexagonal phase than in the 
cubic one. As a first guess, one might attribute this difference to the shorter Eu-Eu distance 
in the hexagonal phase (cIeu-Eu is 3.50-3.74 A in the hexagonal phase and 3.84 A in the cubic 
one). To verify this possibility, we performed test calculations in the hexagonal structure, 
changing the lattice constant so as to match the interatomic Eu-Eu distance of the cubic 
phase. The corresponding exchange coupling did not decrease sufficiently to support this 
view. The rational for the different couplings may probably be found in the different local 
O coordination around Eu ions, as the oxygen orbitals mediate the Eu-Eu interaction: in 
the hexagonal structure there are two O atoms binding two NN Eu ions, while there is only 
one intermediate O atom in the cubic one. An analysis of the partial densities of states 
in the two structures supports this view, with a quite large hybridization of Eu 4/ and O 
orbitals in the hexagonal phase. Interestingly we also observe that the phase with stronger 
exchange coupling has a collinear spin structure, while the structure with weaker interaction 
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is characterized by non-collinear spins. 



VI. CONCLUSION 

In conclusion, we studied the electronic structure of cubic bixbyte EU2O3 and its high 
pressure hexagonal phase by mean of full potential linearized augmented plane wave calcu- 
lations, within the LDA+[7 method. In both phases the filled O 2p and empty majority spin 
Eu 4/ states are separated by a ~ 2.1 eV gap, while minority Eu 4/ states start around 4 
eV above Ep, in agreement with the optical and transport measurements^ 1 ^. 

From the comparison of FM and AFM calculations we infer that the interatomic exchange 
interaction is negligible in cubic EU2O3, and so is its effect on the susceptibility. This is 
consistent with the experimental observation of a constant susceptibility (per mole of Eu) 
of solid solutions of EU2O3 into A2O3 (A=Y, Lu, Sc). Our calculations are in contrast 
with the point of view of Van Vleck and Huang^i, according to which exchange is needed 
to explain the behavior of the magnetic susceptibility. Our results are consistent with 
later calculations^ 3 - which could explain the experimental susceptibility by including matrix 
element in the Van Vleck contribution. In the hexagonal phase we observe the presence 
of a small but sizable AFM exchange interaction between the Eu ion and its six nearest 
neighbors. Therefore, in this phase we could have a non negligible contribution of exchange 
to the magnetic susceptibility, but for now no experimental measurements are available for 
comparison. 
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